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Abstract 

Wireless power transfer (WPT) is a promising new solution to provide convenient and perpetual energy 
supplies to wireless networks. In practice, WPT is implementable by various technologies such as inductive 



I coupling, magnetic resonate coupling, and electromagnetic (EM) radiation, for short-/mid-/long-range applications, 

^ , respectively. In this paper, we consider the EM or radio signal enabled WPT in particular. Since radio signals 

^ I can carry energy as well as information at the same time, a unified study on simultaneous wireless information 



and power transfer (SWIPT) is pursued. Specifically, this paper studies a multiple-input multiple-output (MIMO) 
wireless broadcast system consisting of three nodes, where one receiver harvests energy and another receiver decodes 
information separately from the signals sent by a common transmitter, and all the transmitter and receivers may 
be equipped with multiple antennas. Two scenarios are examined, in which the information receiver and energy 
receiver are separated and see different MIMO channels from the transmitter, or co-located and see the identical 
Q ■ MIMO channel from the transmitter For the case of separated receivers, we derive the optimal transmission strategy 

to achieve different tradeoffs for maximal information rate versus energy transfer, which are characterized by the 

■ boundary of a so-called rate-energy (R-E) region. For the case of co-located receivers, we show an outer bound for 

■ the achievable R-E region due to the potential limitation that practical energy harvesting receivers are not yet able to 



decode information directly. Under this constraint, we investigate two practical designs for the co-located receiver 



■"sj" . case, namely time switching and power splitting, and characterize their achievable R-E regions in comparison to 

. the outer bound. 

o 

Index Terms 



MIMO system, broadcast channel, precoding, wireless power, simultaneous wireless information and power 
transfer (SWIPT), rate-energy ti'adeoff, energy harvesting. 

I. Introduction 

Energy-constrained wireless networks, such as sensor networks, are typically powered by batteries that have 
limited operation time. Although replacing or recharging the batteries can prolong the lifetime of the network to 
a certain extent, it usually incurs high costs and is inconvenient, hazardous (say, in toxic environments), or even 
impossible (e.g., for sensors embedded in building structures or inside human bodies). A more convenient, safer, 
as well as "greener" alternative is thus to harvest energy from the environment, which virtually provides perpetual 
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energy supplies to wireless devices. In addition to other commonly used energy sources such as solar and wind, 
ambient radio-frequency (RF) signals can be a viable new source for energy scavenging. It is worth noting that 
RF-based energy harvesting is typically suitable for low-power applications (e.g., sensor networks), but also can 
be applied for scenarios with more substantial power consumptions if dedicated wireless power transmission is 
implementedy 

On the other hand, since RF signals that carry energy can at the same time be used as a vehicle for transporting 
information, simultaneous wireless information and power transfer (SWIPT) becomes an interesting new area of 
research that attracts increasing attention. Although a unified study on this topic is still in the infancy stage, there 
have been notable results reported in the literature lll]> I0- In Varshney first proposed a capacity-energy function 
to characterize the fundamental tradeoffs in simultaneous information and energy transfer. For the single-antenna 
or SISO (single-input single-output) AWGN (additive white Gaussian noise) channel with amplitude-constrained 
inputs, it was shown in HI that there exist nontrivial tradeoffs in maximizing information rate versus (vs.) power 
transfer by optimizing the input distribution. However, if the average transmit-power constraint is considered instead, 
the above two goals can be shown to be aligned for the SISO AWGN channel with Gaussian input signals, and thus 
there is no nontrivial tradeoff. In [2], Grover and Sahai extended [H to frequency-selective single-antenna AWGN 
channels with the average power constraint, by showing that a non-trivial tradeoff exists in frequency-domain 
power allocation for maximal information vs. energy transfer. 

As a matter of fact, wireless power transfer (WPT) or in short wireless power, which generally refers to the 
transmissions of electrical energy from a power source to one or more electrical loads without any interconnecting 
wires, has been investigated and implemented with a long history. Generally speaking, WPT is carried out using 
either the "near-field" electromagnetic (EM) induction (e.g., inductive coupling, capacitive coupling) for short- 
distance (say, less than a meter) applications such as passive radio-frequency identification (RFID) 131, or the 
"far-field" EM radiation in the form of microwaves or lasers for long-range (up to a few kilometers) applications 
such as the transmissions of energy from orbiting solar power satellites to Earth or spacecrafts |l4]. However, prior 
research on EM radiation based WPT, in particular over the RF band, has been pursued independently from that 
on wireless information transfer (WIT) or radio communication. This is non-surprising since these two lines of 
work in general have very different research goals: WIT is to maximize the information transmission capacity of 
wireless channels subject to channel impairments such as the fading and receiver noise, while WPT is to maximize 
the energy transmission efficiency (defined as the ratio of the energy harvested and stored at the receiver to that 
consumed by the transmitter) over a wireless medium. Nevertheless, it is worth noting that the design objectives 

'interested readers may visit the company website of Powercast at |http://www.powercastco.com/| for more information on recent 
applications of dedicated RF-based power transfer 
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Fig. 1. A wireless network with dual information and energy transfer. 

for WPT and WIT systems can be aligned, since given a transmitter energy budget, maximizing the signal power 
received (for WPT) is also beneficial in maximizing the channel capacity (for WIT) against the receiver noise. 

Hence, in this paper we attempt to pursue a unified study on WIT and WPT for emerging wireless applications 
with such a dual usage. An example of such wireless dual networks is envisaged in Fig. [B where a fixed access 
point (AP) coordinates the two-way communications to/from a set of distributed user terminals (UTs). However, 
unlike the conventional wireless network in which both the AP and UTs draw energy from constant power supplies 
(by e.g. connecting to the grid or a battery), in our model, only the AP is assumed to have a constant power source, 
while all UTs need to replenish energy from the received signals sent by the AP via the far-field RF-based WPT. 
Consequently, the AP needs to coordinate the wireless information and energy transfer to UTs in the downlink, 
in addition to the information transfer from UTs in the uplink. Wireless networks with such a dual information 
and power transfer feature have not yet been studied in the Uterature to our best knowledge, although some of 
their interesting applications have already appeared in, e.g., the body sensor networks ||5] with the out-body local 
processing units (LPUs) powered by battery communicating and at the same time sending wireless power to in-body 
sensors that have no embedded power supplies. However, how to characterize the fundamental information-energy 
transmission tradeoff in such dual networks is still an open problem. 

In this paper, we focus our study on the downlink case with simultaneous WIT and WPT from the AP to UTs. 
In the generic system model depicted in Fig. [H each UT can in general harvest energy and decode information 
at the same time (by e.g. applying the power splitting scheme introduced later in this paper). However, from 
an implementation viewpoint, one particular design whereby each UT operates as either an information receiver 
or an energy receiver at any given time may be desirable, which is referred to as time switching. This scheme 
is practically appealing since state-of-the-art wireless information and energy receivers are typically designed to 
operate separately with very different power sensitivities (e.g., — 50dBm for information receivers vs. — lOdBm 
for energy receivers). As a result, if time switching is employed at each UT jointly with the "near-far" based 
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transmission scheduling at the AP, i.e., UTs that are close to the AP and thus receive high power from the AP are 
scheduled for WET, whereas those that are more distant from the AP and thus receive lower power are scheduled 
for WIT, then SWIPT systems can be efficiently implemented with existing information and energy receivers and 
the additional time-switching device at each receiver. 

For an initial study on SWIPT, this paper considers the simplified scenarios with only one or two active UTs 
in the network at any given time. For the case of two UTs, we assume time switching, i.e., the two UTs take 
tums to receive energy or (independent) information from the AP over different time blocks. As a result, when 
one UT receives information from the AP, the other UT can opportunistically harvest energy from the same signal 
broadcast by the AP, and vice versa. Hence, at each block, one UT operates as an information decoding (ID) 
receiver, and the other UT as an energy harvesting (EH) receiver. We thus refer to this case as separated EH and 
ID receivers. On the other hand, for the case with only one single UT to be active at one time (while all other UTs 
are assumed to be in the off/sleep mode), the active UT needs to harvest energy as well as decode information from 
the same signal sent by the AP, i.e., the same set of receiving antennas are shared by both EH and ID receivers 
residing in the same UT. Thus, this case is referred to as co-located EH and ID receivers. Surprisingly, as we 
will show later in this paper, the optimal information-energy tradeoff for the case of co-located receivers is more 
challenging to characterize than that for the case of separated receivers, due to a potential Umitation that practical 
EH receiver circuits are not yet able to decode the information directly and vice versa. Note that similar to the 
case of separated receives, time switching can also be appUed in the case of co-located receivers to orthogonahze 
the information and energy transmissions at each receiving antenna; however, this scheme is in general suboptimal 
for the achievable rate-energy tradeoffs in the case of co-located receivers, as will be shown later in this paper. 

Some further assumptions are made in this paper for the purpose of exposition. Firstly, this paper considers a 
quasi-static fading enviroimient where the wireless channel between the AP and each UT is assumed to be constant 
over a sufficiently long period of time during which the number of transmitted symbols can be approximately 
regarded as being infinitely large. Under this assumption, we further assume that it is feasible for each UT to 
estimate the downhnk channel from the AP and then send it back to the AP via the uphnk, since the time overhead 
for such chaimel estimation and feedback is a neghgible portion of the total transmission time due to quasi-static 
fading. We will address the more general case of fading channels with imperfect/partial channel knowledge at the 
transmitter in our future work. Secondly, we assume that the system under our study typically operates at the high 
signal-to-noise ratio (SNR) regime for the ID receiver in the case of co-located receivers. This is to be compatible 
with the high-power operating requirement for the EH receiver of practical interest as previously mentioned. Thirdly, 
without loss of generaUty, we assume a multi-antenna or MIMO (multiple-input multiple-output) system, in which 
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Fig. 2. A MIMO broadcast system for simultaneous wireless information and power transfer. 

the AP is equipped with multiple antennas, and each UT is equipped with one or more antennas, for enabling 
both the high-performance wireless energy and information transmissions (as it is well known that for WIT only, 
MIMO systems can achieve folded array/capacity gains over SISO systems by spatial beamforming/multiplexing 
0). 

Under the above assumptions, a three-node MIMO broadcast system is considered in this paper, as shown in 
Fig. m wherein the EH and ID receivers harvest energy and decode information separately from the signal sent 
by a common transmitter. Note that this system model refers to the case of separated EH and ID receivers in 
general, but includes the co-located receivers as a special case when the MIMO channels from the transmitter to 
both receivers become identical. Assuming this model, the main results of this paper are summarized as follows: 

• For the case of separated EH and ID receivers, we design the optimal transmission strategy to achieve different 
tradeoffs between maximal information rate vs. energy transfer, which are characterized by the boundary 
of a so-called rate-energy (R-E) region. We derive a semi-closed-form expression for the optimal transmit 
covariance matrix (for the joint precoding and power allocation) to achieve different rate-energy pairs on the 
boundary of the R-E region. Note that the R-E region is a multiuser extension of the single-user capacity- 
energy function in [1]. Also note that the multi-antenna broadcast channel (BC) has been investigated in e.g. 
|[7l- llT2l for information transfer solely by unicasting or multicasting. However, MIMO-BC for SWIPT as 
considered in this paper is new and has not yet been studied by any prior work. 

• For the case of co-located EH and ID receivers, we show that the proposed solution for the case of separated 
receivers is also applicable with the identical MIMO channel from the transmitter to both ID and EH receivers. 
Furthermore, we consider a potential practical constraint that EH receiver circuits cannot directly decode the 
information (i.e., any information embedded in received signals sent to the EH receiver is lost during the EH 
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Fig. 3. Two practical designs for the co-located energy and information receivers, which are applied for each receiving antenna. 



process). Under this constraint, we show that the R-E region with the optimal transmit covariance (obtained 
without such a constraint) in general only serves as a performance outer bound for the co-located receiver 
case. 

• Hence, we investigate two practical receiver designs, namely time switching and power splitting, for the case 
of co-located receivers. As shown in Fig. |3] for time switching, each receiving antenna periodically switches 
between the EH receiver and ID receiver, whereas for power splitting, the received signal at each antenna is 
split into two separate signal streams with different power levels, one sent to the EH receiver and the other 
to the ID receiver. Note that time switching has also been proposed in [il4il for the SISO AWGN channel. 
Furthermore, note that the antenna switching scheme whereby the receiving antennas are divided into two 
groups with one group switched to information decoding and the other group to energy harvesting can be 
regarded as a special case of power splitting with only binary splitting power ratios at each receiving antenna. 
For these practical receiver designs, we derive their achievable R-E regions as compared to the R-E region 
outer bound, and characterize the conditions under which their performance gaps can be closed. For example, 
we show that the power splitting scheme approaches the tradeoff upper bound asymptotically when the RF- 
band antenna noise at the receiver becomes more dominant over the baseband processing noise (more details 
are given in Section IV-C). 

The rest of this paper is organized as follows: Section |Il] presents the system model, characterizes the rate-energy 
region, and formulates the problem for finding the optimal transmit covariance matrix. Section |IlI] presents the 
optimal transmit covariance solution for the case of separated receivers. Section JV] extends the solution to the case 
of co-located receivers to obtain a performance upper bound, proposes practical receiver designs, and analyzes 
their performance limits as compared to the performance upper bound. Finally, Section |V] concludes the paper and 
provides some promising directions for future work. 

Notation: For a square matrix S, tr(S'), l^], S^^, and 5^ denote its trace, determinant, inverse, and square- 
root, respectively, while S ^ and S ^ mean that S is positive semi-definite and positive definite, respectively. 
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For an arbitrary-size matrix M, and M'^ denote the conjugate transpose and transpose of M, respectively. 
diag(xi, . . . ,xm) denotes an M x M diagonal matrix with xi,. . . ,xm being the diagonal elements. / and 
denote an identity matrix and an all-zero vector, respectively, with appropriate dimensions. E[-] denotes the statistical 
expectation. The distribution of a circularly symmetric complex Gaussian (CSCG) random vector with mean x and 
covariance matrix S is denoted by CJ\f{x, XI), and ~ stands for "distributed as". C^^^ denotes the space of x x y 
matrices with complex entries. ||z|| is the Euclidean norm of a complex vector z, and \z\ is the absolute value of 
a complex scalar z. max(x,y) and mm{x,y) denote the maximum and minimum between two real numbers, x 
and y, respectively, and (x)^ = max(j;,0). All the log(-) functions have base-2 by default. 

II. System Model And Problem Formulation 

As shown in Fig. [21 this paper considers a wireless broadcast system consisting of one transmitter, one EH 
receiver, and one ID receiver. It is assumed that the transmitter is equipped with M > 1 transmitting antennas, and 
the EH receiver and the ID receiver are equipped with A'^eh > 1 and Nid > 1 receiving antennas, respectively. In 
addition, it is assumed that the transmitter and both receivers operate over the same frequency band. Assuming a 
narrow-band transmission over quasi-static fading channels, the baseband equivalent channels from the transmitter 
to the EH receiver and ID receiver can be modeled by matrices G € C^'^'^^^^ and H E qNidxM ^ respectively. 
It is assumed that at each fading state, G and H are both known at the transmitter, and separately known at the 
corresponding receiver. Note that for the case of co-located EH and ID receivers, G is identical to H and thus 
Ner = Nib. 

It is worth noting that the EH receiver does not need to convert the received signal from the RF band to the 
baseband in order to harvest the carried energy. Nevertheless, thanks to the law of energy conservation, it can be 
assumed that the total harvested RF-band power (energy normalized by the baseband symbol period), denoted by 
Q, from all receiving antennas at the EH receiver is proportional to that of the received baseband signal, i.e., 

Q = CmGx{n)f] (1) 

where C is a constant that accounts for the loss in the energy transducer for converting the harvested energy to 
electrical energy to be stored; for the convenience of analysis, it is assumed that ^ = 1 in this paper unless stated 
otherwise. We use x{n) G c^^xi to denote the baseband signal broadcast by the transmitter at the nth symbol 
interval, which is assumed to be random over n, without loss of generality. The expectation in ([T]) is thus used to 
compute the average power harvested by the EH receiver at each fading state. Note that for simplicity, we assumed 
in ([Til that the harvested energy due to the background noise at the EH receiver is negligible and thus can be 
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ignoredy 

On the other hand, the baseband transmission from the transmitter to the ID receiver can be modeled by 

y{n) = Hx{n) + z{n) (2) 

where y{n) G c^mxi (jguotes the received signal at the nth symbol interval, and z{n) G c^mxi denotes the 
receiver noise vector. It is assumed that z(n)'s are independent over n and z{n) ^ CJ\f{0, 1). Under the assumption 
that x{n) is random over n, we use S = ¥,[x{n)x^ (n)] to denote the covariance matrix of x{n). In addition, 
we assume that there is an average power constraint at the transmitter across all transmitting antennas denoted 
by E[||a;(n)|p] = tr{S) < P. In the following, we examine the optimal transmit covariance S to maximize the 
transported energy efficiency and information rate to the EH and ID receivers, respectively. 

Consider first the MIMO link from the transmitter to the EH receiver when the ID receiver is not present. In 
this case, the design objective for S is to maximize the power Q received at the EH receiver. Since from (H) it 
follows that Q = tr{GSG^) with C = 1, the aforementioned design problem can be formulated as 

(PI) max Q := tr {GSG^) 

s 

s.t. tr(5)<P,5^0. 

Let Ti = min(M, A^eh) and the (reduced) singular value decomposition (SVD) of G be denoted by G = 
UgTq'^Vq, where Uq G C^ehxTi ^ C*^^^\ each of which consists of orthogonal columns with 

unit norm, and Tg = diag((jri, . . . , (^^J with gi > §2 > ■ ■ ■ > QTi > 0. Furthermore, let i^i denote the first 
column of Vq- Then, we have the following proposition. 

Proposition 2.1: The optimal solution to (PI) is S^eh = Pviv^ . 

Proof: See Appendix lAl ■ 

Given S = Seu, it follows that the maximum harvested power at the EH receiver is given by Qmax = 9iP- It is 
worth noting that since Seu is a rank-one matrix, the maximum harvested power is achieved by beamforming at 
the transmitter, which aligns with the strongest eigenmode of the matrix G^G, i.e., the transmitted signal can be 
written as x{n) = \/Pvis{n), where s{n) is an arbitrary random signal over n with zero mean and unit variance, 
and vi is the transmit beamforming vector. For convenience, we name the above transmit beamforming scheme 
to maximize the efficiency of WPT as "energy beamforming". 

Next, consider the MIMO link from the transmitter to the ID receiver without the presence of any EH receiver. 

Assuming the optimal Gaussian codebook at the transmitter, i.e., x{n) ~ CA/'(0, S), the transmit covariance S to 

^The results of this paper are readily extendible to study the impacts of non-negligible background noise and/or co-channel interference 
on the SWIPT system performance. 
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maximize the transmission rate over this MIMO channel can be obtained by solving the following problem |fT3l|: 



The optimal solution to the above problem is known to have the following form [13 1: Sm = Vh-^^H' where 
G C^^^"^^ is obtained from the (reduced) SVD of H expressed hy H = jy^r^Vg, with T2 = min(M, iVio), 
Uh e C^i^x'T^ = diag(/ii,...,/iTj, /ii > /i2 > ... > /it. > 0, and A = diag(pi, . . . ,ptJ with the 
diagonal elements obtained from the standard "water-filling (WF)" power allocation solution llT3l : 



with u being the so-called (constant) water-level that makes J2iliPi — The corresponding maximum trans- 
mission rate is then given by i?max = Yll=i log(l + hiPi). The maximum rate is achieved in general by spatial 
multiplexing IS over up to T2 spatially decoupled AWGN channels, together with the Gaussian codebook, i.e., the 



Remark 2.1: It is worth noting that in Problem (PI), it is assumed that the transmitter sends to the EH receiver 
continuously. Now suppose that the transmitter only transmits a fraction of the total time denoted by a with 
< a < 1. Furthermore, assume that the transmit power level can be adjusted flexibly provided that the consumed 
average power is bounded by P, i.e., a ■ tr(S) + (1 — a) . < P or tr(S) < P/a. In this case, it can be easily 
shown that the transmit covariance S = {P / a)vivf also achieves the maximum harvested power Qmax = QiP for 
any < a < 1, which suggests that the maximum power delivered is independent of transmission time. However, 
unlike the case of maximum power transfer, the maximum information rate reliably transmitted to the ID receiver 
requires that the transmitter send signals continuously, i.e., a = 1, as assumed in Problem (P2). This can be 
easily verified by observing that for any < a < 1 and 5^0, a log \I + H{S/a)H^\ < log \I + HSH^\ 
where the equality holds only when a = 1, since R is a nonlinear concave function of S. Thus, to maximize 
both power and rate transfer at the same time, the transmitter should broadcast to the EH and ID receivers all the 
time. Furthermore, note that the assumed Gaussian distribution for transmitted signals is necessary for achieving 
the maximum rate transfer, but not necessary for the maximum power transfer. In fact, for any arbitrary complex 
number c that satisfies |c| = 1, even a deterministic transmitted signal x{n) = ^/Pvic,\/n, achieves the maximum 
transferred power Qmax in Problem (PI). However, to maximize simultaneous power and information transfer with 
the same transmitted signal, the Gaussian input distribution is sufficient as well as necessary. 

Now, consider the case where both the EH and ID receivers are present. From the above results, it is seen that the 



(P2) max R:=log\I + HSH 

s 

s.t. tr{S)<P,ShO. 




(3) 



transmitted signal can be expressed as x{n) = Vi/A^/^s(n), where s{n) is a Gaussian random vector ~ CM{0, 1), 
V H and A^/^ denote the precoding matrix and the (diagonal) power allocation matrix, respectively. 
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Fig. 4. Rate-energy tradeoff for a MIMO broadcast system with separated EH and ID receivers, and M = A'^eh ~ Nio — 4. 



optimal transmission strategies for maximal power transfer and information transfer are in general different, which 
are energy beamforming and information spatial multiplexing, respectively. It thus motivates our investigation of 
the following question: What is the optimal broadcasting strategy for simultaneous wireless power and information 
transfer? To answer this question, we propose to use the Rate-Energy (R-E) region (defined below) to characterize 
all the achievable rate (in bits/sec/Hz or bps for information transfer) and energy (in joule/sec or watt for power 
transfer) pairs under a given transmit power constraint. Without loss of generality, assuming that the transmitter 
sends Gaussian signals continuously (cf. Remark 12.11 ). the R-E region is defined as 



Cr_e(P) = : R < log\I + HSH^\,Q < tr(G5G^),tr(5) < P,Sh o|. 



(4) 



In Fig. m an example of the above defined R-E region (see Section |lll] for the algorithm to compute the boundary 
of this region) is shown for a practical MIMO broadcast system with separated EH and ID receivers (i.e., G ^ H). 
It is assumed that M = A'^eh = -^ID = 4. The transmitter power is assumed to be P = 1 watt(W) or 30dBm. The 
distances from the transmitter to the EH and ID receivers are assumed to be 1 meter and 10 meters, respectively; 
thus, we can exploit the near-far based energy and information transmission scheduling, which may correspond 
to, e.g., a dedicated energy transfer system (to "near" users) with opportunistic information transmission (to "far" 
users), or vice versa. Assuming a carrier frequency of /c = 900MHz and the power pathloss exponent to be 4, 
the distance-dependent signal attenuation from the AP to EH/ID receiver can be estimated as 40dB and 80dB, 
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respectively. Accordingly, the average signal power at the EH/ID receiver is thus 30dBm— 40dB= — lOdBm and 
30dBm— 80dB= — 50dBm, respectively. It is further assumed that in addition to signal pathloss, Rayleigh fading is 
present, as such each element of channel matrices G and H is independently drawn from the CSCG distribution with 
zero mean and variance — lOdBm (for EH receiver) and — 50dBm for (for ID receiver), respectively (to be consistent 
with the signal pathloss previously assumed). Furthermore, the bandwidth of the transmitted signal is assumed to 
be lOMHz, while the receiver noise is assumed to be white Gaussian with power spectral density — 140dBm/Hz 
(which is dominated by the receiver processing noise rather than the background thermal noise) or average power 
— 70dBm over the bandwidth of lOMHz. As a result, considering all of transmit power, signal attenuation, fading 
and receiver noise, the per-antenna average SNR at the ID receiver is equal to 30 — 80 — (—70) = 20dB, which 
corresponds to P = 100 in the equivalent signal model for the ID receiver given in with unit-norm noise. 
In addition, we assume that for the EH receiver, the energy conversion efficiency is C =50%. Considering this 
together with transmit power and signal attenuation, the average per-antenna signal power at the EH receiver is 
thus 0.5x(30dBm-40dB) = 50/xW. 

From Fig. IH it is observed that with energy beamforming, the maximum harvested energy rate for the EH 
receiver is around Qmax = 0.57mW, while with spatial multiplexing, the maximum information rate for the ID 
receiver is around i?max = 225Mbps. It is easy to identify two boundary points of this R-E region denoted by 
{Reh, Qmax) and (i?max) Qid), respectively. For the former boundary point, the transmit co variance is Seh, which 
corresponds to transmit beamforming and achieves the maximum transferred power Qmax to the EH receiver, while 
the resulting information rate for the ID receiver is given by Rer = log(l + pP). On the other hand, for 

the latter boundary point, the transmit co variance is Sm, which corresponds to transmit spatial multiplexing and 
achieves the maximum information rate transferred to the ID receiver i?max> while the resulting power transferred 
to the EH receiver is given by Qid = tr(GS'iDG^). 

Since the optimal tradeoff between the maximum energy and information transfer rates is characterized by the 
boundary of the R-E region, it is important to characterize all the boundary rate-power pairs of Cr_e(^) for any 
P > 0. From Fig. |4] it is easy to observe that if i? < Rer, the maximum harvested power Qmax is achievable with 
the same transmit covariance that achieves the rate-power pair (Peh, Qmax); similarly, the maximum information 
rate Pmax is achievable provided that Q < Qid. Thus, the remaining boundary of Cr_e(P) yet to be characterized 
is over the intervals: Rer < R < Rma.x, Qid < Q < Qmax- We thus consider the following optimization problem: 

(P3) max log\l + HSH^l 

s 

s.t. tr [GSG^) > Q, tr(5) < P, StO- 
Note that if Q takes values from Qid < Q < Qmax> the corresponding optimal rate solutions of the above problems 
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are the boundary rate points of the R-E region over Reh < R < Rma.x- Notice that the transmit covariance solutions 
to the above problems in general yield larger rate-power pairs than those by simply "time-sharing" the optimal 
transmit covariance matrices Seu and Sib for EH and ID receivers separately (see the dashed line in Fig. ©El 

Problem (P3) is a convex optimization problem, since its objective function is concave over S and its constraints 
specify a convex set of S. Note that (P3) resembles a similar problem formulated in ||T5l . |[T6l (see also fVJ] and 
references therein) under the cognitive radio (CR) setup, where the rate of a secondary MIMO Unk is maximized 
subject to a set of so-called interference power constraints to protect the co-channel primary receivers. However, 
there is a key difference between (P3) and the problem in ||T61 : the harvested power constraint in (P3) has the 
reversed inequality of that of the interference power constraint in fl^, since in our case it is desirable for the 
EH receiver to harvest more power from the transmitter, as opposed to that in 11161 the interference power at the 
primary receiver should be minimized. As such, it is not immediately clear whether the solution in |[T6l can be 
directly applied for solving (P3) with the reversed power inequality. In the following, we will examine the solutions 
to Problem (P3) for the two cases with arbitrary G and H (the case of separated receivers) and G = H (the case 
of co-located receivers), respectively. 

III. Separated Receivers 

Consider the case where the EH receiver and ID receiver are spatially separated and thus in general have different 
channels from the transmitter. In this section, we first solve Problem (P3) with arbitrary G and H and derive a 
semi-closed-form expression for the optimal transmit covariance. Then, we examine the optimal solution for the 
special case of MISO channels from the transmitter to ID and/or EH receivers. 

Since Problem (P3) is convex and satisfies the Slater's condition ifTSl . it has a zero duality gap and thus can 
be solved using the Lagrange duality methodo Thus, we introduce two non-negative dual variables, A and /x, 
associated with the harvested power constraint and transmit power constraint in (P3), respectively. The optimal 
solution to Problem (P3) is then given by the following theorem in terms of A* and /x*, which are the optimal 
dual solutions of Problem (P3) (see Appendix |B] for details). Note that for Problem (P3), given any pair of Q 
(QiD < Q < Qmax) and P > 0, there exists one unique pair of A* > and /i* > 0. 

Theorem 3.1: The optimal solution to Problem (P3) has the following form: 

S* = A-^l'^VkV^ A-^l'^ (5) 

^By time-sharing, we mean that the AP transmits simultaneously to both EH and ID receivers with the energy-maximizing transmit 
covariance Seh (i.e. energy beamforming) for /3 portion of each block time, and the information-rate-maximizing transmit covariance Syd 
(i.e. spatial multiplexing) for the remaining 1 — /? portion of each block time, with < /? < 1. 

''it is worth noting that Problem (P3) is convex and thus can be solved efficiently by the interior point method III8I : in this paper, we 
apply the Lagrange duality method for this problem mainly to reveal the optimal precoder structure. 
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where A = fi*I - X*G^G, V £ C^^^T2 obtained from the (reduced) SVD of the matrix iJA"^/^ gjygn by 
HA-^/'^ = ilt^^'^V^, with f = diag(^i, . . . ,/itJ, /ii > /i2 > • • • > /it. > 0, and A = diag(pi, . . . ,prj, 

with Pi = {l-l/hi)+,i = l,...,T2. 



Note that this theorem requires that A = - X*G G >~ 0, implying that fi* > X*gi (recall that gi is the 
largest eigenvalue of matrix G^G), which is not present for a similar result in [17] under the CR setup with 
the reversed interference power constraint. One algorithm that can be used to solve (P3) is provided in Table 
|I] of Appendix |Bj From Theorem 13.11 the maximum transmission rate for Problem (P3) can be shown to be 
R* = log |/ + HS*H^\ = YldLi log(l + hiPi), for which the proof is omitted here for brevity. 

Next, we examine the optimal solution to Problem (P3) for the special case where the ID receiver has one single 
antenna, i.e., Nm = 1, and thus the MIMO channel H reduces to a row vector with h G C^^^^. Suppose 
that the EH receiver is still equipped with A'^eh > 1 antennas, and thus the MIMO channel G remains unchanged. 
From Theorem 13.11 we obtain the following corollary. 

Corollary 3.1: In the case of MISO channel from the transmitter to ID receiver, i.e., H = , the optimal 
solution to Problem (P3) reduces to the following form: 



where A = fJ,*I — X*G G, with A* and fi* denoting the optimal dual solutions of Problem (P3). Correspondingly, 
the optimal value of (P3) is R* = (2 log (\\A-^/'^h\\)Y . 



From it is observed that the optimal transmit covariance is a rank-one matrix, from which it follows that 
beamforming is the optimal transmission strategy in this case, where the transmit beamforming vector should be 
aligned with the vector A^h. Moreover, consider the case where both channels from the transmitter to ID/EH 
receivers are MISO, i.e., H = , and G = with g G <[^Mxi^ From Corollary 13.11 it follows immediately 
that the optimal covariance solution to Problem (P3) is still beamforming. In the following theorem, we show a 
closed-form solution of the optimal beamforming vector at the transmitter for this special case, which differs from 
the semi-closed-form solution © that was expressed in terms of dual variables. 

Theorem 3.2: In the case of MISO channels from transmitter to both ID and EH receivers, i.e., H = h^, and 
G = g^ , the optimal solution to Problem (P3) can be expressed as S* = Pvv^ , where the beamforming vector 
V has a unit-norm and is given by 



Proof: See Appendix IbI 




(6) 





Proof: See Appendix ICl 
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Fig. 5. Rate-energy tradeoff for a MISO broadcast system with coixelated MISO channels to the (separated) EH and ID receivers. 

where h = h/\\h\\, g = g/\\g\\, hg± = hg±/\\hg±\\ with hg± = h-{g"h)g, and Ugh = g^h with Zagh G [0, 27r) 
denoting the phase of complex number agh- Correspondingly, the optimal value of (P3) is given by 



R* 



log(l + \\hfP) 
loK (1 + 



< Q < \g^h\'^P 



g^hl'^P <Q< P\\g\\'^. 



(8) 



Proof: The proof is similar to that of Theorem 2 in [16|, and is thus omitted for brevity. ■ 
It is worth noting that in (O, if Q < \g^h\'^P, the optimal transmit beamforming vector is based on the 
principle of maximal-ratio-combining (MRC) with respect to the MISO channel from the transmitter to the 
ID receiver, and in this case, the harvested power constraint in Problem (P3) is indeed not active; however, when 
Q > \g^h\'^P, the optimal beamforming vector is a linear combination of the two vectors g and hg±, and the 
combining coefficients are designed such that the harvested power constraint is satisfied with equality. 

In Fig. 121 we show the achievable R-E regions for the case of MISO channels from the transmitter to both 
EH and ID receivers. We set P = 10. For the purpose of exposition, it is assumed that \\h\\ = \\g\\ = 1 and 
lo^g/iP = P» with < p < 1 denoting the correlation between the two unit-norm vectors h and g. This channel 
setup may correspond to the practical scenario where the EH and ID receivers are equipped at a single device 
(but still physically separated), and as a result their respective MISO channels from the transmitter have the same 
power gain but are spatially correlated due to the insufficient spacing between two separate receiving antennas. 
From Theorem 13.21 the R-E regions for the three cases of p = 0.1,0.5, and 0.9 are obtained, as shown in Fig. |5] 
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Interestingly, it is observed that increasing p enlarges the achievable R-E region, which indicates that the antenna 
correlation between the EH and ID receivers can be a beneficial factor for simultaneous information and power 
transfer. Note that in this figure, we express energy and rate in terms of energy unit and bits/channel use, respectively, 
since their practical values can be obtained by appropriate scaling based on the realistic system parameters as for 



In this section, we address the case where the EH and ID receivers are co-located, and thus possess the same 
channel from the transmitter, i.e., G = H and thus A'^eh = ^ID — N. We first examine the optimal solution of 
Problem (P3) for this case, from which we obtain an outer bound for the achievable rate-power pairs in the R-E 
region. Then, we propose two practical receiver designs, namely time switching and power splitting, derive their 
optimal transmission strategies to maximize the achievable rate -power pairs, and finally compare the results to the 
R-E region outer bound. 

A. Performance Outer Bound 

Consider Problem (P3) with G = H. Recall that the (reduced) SVD of H is given by ii" = U h^^h^Vh, with 
Vh = diag(ft,i, . . . , hx^), hi > h2 > ■ ■ ■ > h^^ > 0, and T2 = min(M, N). From Theorem 13.11 we obtain the 
following corollary. 

Corollary 4.1: In the case of co-located EH and ID receivers with G = H, the optimal solution to Problem 
(P3) has the form of S* = Vh'^V^, where Xl = diag(pi, . . . jPTa) with the diagonal elements obtained from 
the following modified WF power allocation: 



with A* and n* denoting the optimal dual solutions of Problem (P3), fi* > X*hi. The corresponding maximum 
transmission rate is R* = log(l + KPi)- 



The algorithm in Table U for solving Problem (P3) with arbitrary G and H can be simplified to solve the special 
case with G = H. Corollary 14. 1 1 reveals that for Problem (P3) in the case of G = H, the optimal transmission 
strategy is in general still spatial multiplexing over the eigenmodes of the MIMO channel H as for Problem (P2), 
while the optimal tradeoffs between information transfer and power transfer are achieved by varying the power 
levels allocated into different eigenmodes, as shown in dill. It is interesting to observe that the power allocation 
in ^ reduces to the conventional WF solution in ^ with a constant water-level when A* = 0, i.e., the harvested 
power constraint in Problem (P3) is inactive with the optimal power allocation. However, when A* > and thus 



Fig. a 



IV. Co-Located Receivers 




(9) 



Proof: See Appendix iDl 
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the harvested power constraint is active corresponding to the Pareto-optimal regime of our interest, the power 
allocation in ([91) is observed to have a non-decreasing water-level as /ij's increase. Note that this modified WF 
policy has also been shown in [2] for power allocation in frequency-selective AWGN channels. 

Using Corollary 14. II we can characterize all the boundary points of the R-E region Cji-e{P) defined in Q for the 
case of co-located receivers with G = H. For example, if the total transmit power is allocated to the channel with 
the largest gain hi, i.e., pi = P and pi = 0,i = 2, . . . ,T2, the maximum harvested power Qmax = Phi is achieved 
by transmit beamforming. On the other hand, if transmit spatial multiplexing is applied with the conventional WF 
power allocation given in (|9]l with A* = 0, the corresponding R* becomes the maximum transmission rate, i?max- 
However, unlike the case of separated EH and ID receivers in which the entire boundary of Cr_e(^) is achievable, 
in the case of co-located receivers, except the two boundary rate-power pairs (-Rmaxi 0) and (0, Qma.x), all the other 
boundary pairs of Cr„e(^) may not be achievable in practice. Note that these boundary points are achievable 
if and only if (iff) the following premise is true: the power of the received signal across all antennas is totally 
harvested, and at the same time the carried information with a transmission rate up to the MIMO channel capacity 
(for a given transmit covariance) is decodable. However, existing EH circuits are not yet able to directly decode 
the information carried in the RF-band signal, even for the SISO channel case; as a result, how to achieve the 
remaining boundary rate-power pairs of Cr_e(^) in the MIMO case with the co-located EH and ID receiver 
remains an interesting open problem. Therefore, in the case of co-located receivers, the boundary of Cr_e(^) 
given by Corollary 14.11 in general only serves as an outer bound for the achievable rate -power pairs with practical 
receiver designs, as will be investigated in the following subsections. 

B. Time Switching 

First, as shown in Fig. Ha), we consider the time switching (TS) scheme, with which each transmission block 
is divided into two orthogonal time slots, one for transferring power and the other for transmitting data. The co- 
located EH and ID receiver switches its operations periodically between harvesting energy and decoding information 
between the two time slots. It is assumed that time synchronization has been perfectly established between the 
transmitter and the receiver, and thus the receiver can synchronize its function switching with the transmitter. With 
orthogonal transmissions, the transmitted signals for the EH receiver and ID receiver can be designed separately, 
but subject to a total transmit power constraint. Let a with < a < 1 denote the percentage of transmission time 
allocated to the EH time slot. We then consider the following two types of power constraints at the transmitter: 

• Fixed power constraint: The transmitted signals to the ID and EH receivers have the same fixed power 
constraint given by tr(S^i) < P, and tr(S'2) < P, where S^i and ^2 denote the transmit covariance matrices 
for the ID and EH transmission time slots, respectively. 
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• Flexible power constraint: The transmitted signals to the ID and EH receivers can have different power 
constraints provided that their average consumed power is below P, i.e., (1 — a)tr(S'i) + atr(S2) < P- 

Note that the TS scheme under the fixed power constraint has been considered in |14| for the single-antenna 
AWGN channel. The achievable R-E regions for the TS scheme with the fixed (referred to as TSi) vs. flexible 
(referred to as TS2) power constraints are then given as follows: 



Cl%^P)= U <(!-«) log 

0<a<l ^ 



\I + HSiH"l 



Q < atr{HS2H"),tr{Si) < P, tr(52) < P \ (10) 



Cr!e(^)= U {iR,Q)-R<i^-(^nog\i + HSiH''\, 

0<a<l ^ 

Q < atr{HS2H"), (1 - a)tr(Si) + atr(52) < p|. (11) 

It is worth noting that C'^%{P) C C^^"p.(P) must be true since any pair of 5i ^ and ^2 ^ that satisfy the 
fixed power constraint will satisfy the flexible power constraint, but not vice versa. The optimal transmit covariance 
matrices Si and ^2 to achieve the boundary of Cj^^g(P) with the fixed power constraint are given in Section HI] 
(assuming G = H). In fact, the boundary of C^-^{P) is simply a straight line connecting the two points (i?max> 0) 
and (0, Qmax) (cf. Fig. IT]) by sweeping a from to 1. 

Similarly, for the case of flexible power constraint, the transmit covariance solutions for Si and ^2 to achieve 
any boundary point of C^J^{P) can be shown to have the same set of eigenvectors as those given in Section HIl 
(assuming G = H), respectively; however, the corresponding time allocation for a and power allocation for Si 
and S2 remain unknown. We thus have the following proposition. 

Proposition 4.1: In the case of flexible power constraint, except the two points (PmaxjO) and (0,(5max)> all 
Other boundary points of the region C-^J^{P) are achieved as a — )■ 0; accordingly, C-^J^{P) can be simplified as 

Cl%{P) = \{R,Q):R<\og\I + HSiH''ltx{Si) < (P - Q/Zii), 5i ^ o| . (12) 

Proof: See Appendix IE] ■ 
The corresponding optimal power allocation for Si and S2 can be easily obtained given ([T2l ) and are thus omitted 
for brevity. Proposition 14. 1 1 suggests that to achieve any boundary point {R,Q) of Cj^^g(P) with R < i?max and 
Q < QmsLx, the portion of transmission time a allocated to power transfer in each block should asymptotically 
go to zero when n — )• 00, where n denotes the number of transmitted symbols in each block. For example, by 
allocating O(logn) symbols per block for power transfer and the remaining symbols for information transmission 
yields a = log n/n — )• as n — )• 00, which satisfies the optimality condition given in Proposition 14. 1 1 
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Fig. 6. Receiver operations with/without a power splitter (the energy harvested due to the receiver noise is ignored for EH receiver). 

It is worth noting that the boundary of C^^{P) in the flexible power constraint case is achieved under the 
assumption that the transmitter and receiver can both operate in the regime of infinite power in the EH time slot 
due to a — > 0, which cannot be implemented with practical power amplifiers. Hence, a more feasible region for 
Cj^^g(P) is obtained by adding peajf] transmit power constraints in (ITTI i as tr(S^i) < Ppcak and tr(S^2) < Ppcak> 
with Ppcak ^ P- Similar to Proposition 14.11 it can be shown that the boundary of the achievable R-E region in 
this case, denoted by C^^{P, Ppcak)> is achieved by a = Q/(/iiPpcak)- Note that we can equivalently denote the 
achievable R-E region Cr!e(P) defined in ([n) or (O without any peak power constraint as C^^g(P, oo). 

C. Power Splitting 

Next, we propose an alternative receiver design called power splitting (PS), whereby the power and information 
transfer to the co-located EH and ID receivers are simultaneously achieved via a set of power splitting devices, 
one device for each receiving antenna, as shown in Fig. Sb). In order to gain more insight into the PS scheme, we 
consider first the simple case of a single-antenna AWGN channel with co-located ID and EH receivers, which is 
shown in Fig. Sa). For the ease of comparison, the case of solely information transfer with one single ID receiver 
is also shown in Fig. I^b). 

The receiver operations in Fig. |6la) are explained as follows: The received signal from the antenna is first 
corrupted by a Gaussian noise denoted by nA{t) at the RF-band, which is assumed to have zero mean and 
equivalent baseband power a\. The RF-band signal is then fed into a power splitter, which is assumed to be 
perfect without any noise induced. After the power splitter, the portion of signal power split to the EH receiver is 

'Note that the peak power constraint in this context is different from the signal amplitude constraint considered in IT], 1141 . 



19 

denoted by p, and that to the ID receiver by 1 — p. The signal split to the ID receiver then goes through a sequence 
of standard operations (see, e.g., iflQl ) to be converted from the RF band to baseband. During this process, the 
signal is additionally corrupted by another noise np{t), which is independent of nA{t) and assumed to be Gaussian 
and have zero mean and variance Up. To be consistent to the case with solely the ID receiver, it is reasonable to 
assume that the antenna noise n^(t) and processing noise np{t) have the same distributions in both Figs. |6la) and 
[6tb). It is further assumed that cj^ + dp = 1 to be consistent with the system model introduced in Section Ull 

For this simple SISO AWGN channel, we denote the transmit power constraint by P and the channel power 
gain by h. It is then easy to compute the R-E region outer bound Cr-eI^') for this channel with co-located ID 
and EH receivers, which is simply a box specified by three vertices (0, (5max)> (^max> 0) and (iimaxi Qmax)> with 
Qmax = Ph and i^max = log(l + Ph). Interestingly, we will show next that under certain conditions, the PS 
scheme can in fact achieve all the rate-energy pairs in this R-E region outer bound; without loss of generality, it 
suffices to show that the vertex point (iimax, Qmax) is achievable. 

With reference to Fig. IHa), we discuss the PS scheme in the following three regimes with different values of 
antenna and processing noise power. 

• <t\ <^ ap (Case I): In this ideal case with perfect receiving antenna, the antenna noise can be ignored and 
thus we have a\ = and ap = 1. Accordingly, it is easy to show that the SNR, denoted by r, at the ID 
receiver in Fig. |6la) is (1 — p)Ph. The achievable R-E region in this case is then given by C^^'g(P) = 
Uo<p<i{(^' ^3) • ^ — los(l + (1 ~ p)Ph),Q < pPh}. This region can be shown to coincide with the R-E 
region for the TS scheme with the flexible power constraint given by (fT2l) for the SISO case. 

• < o"^ < 1 (Case II): This is the most practically valid case. Since dp = 1 — a\, we can show that r in this 
case is given by r = (1 — p)Ph/ ((1 — p)(y\ + dp) = (1 — p)Ph/ (1 — pa\). Accordingly, the achievable R-E 
region in this case is given by C^^'^^(P) = Uo<p<ii(-^' Q) '■ P ^ log(l + t)-,Q < pPh}. It is easy to show 
that C^^'g (P) enlarges strictly as a\ increases from to 1. 

• d^ ^ dp (Case III): In this ideal case with perfect RF-to-baseband signal conversion, the processing noise can 
be ignored and thus we have dp = and ct\ = 1. In this case, the SNR for the ID receiver is given by r = Ph, 
regardless of the value of p. Thus, to maximize the power transfer, ideally we should set p — 1, i.e., splitting 
infinitesimally small power to the ID receiver since both the signal and antenna noise are scaled identically by 
the power splitter and there is no additional processing noise induced after the power splitting. With p = I, 
the achievable R-E region in this case is given by C^_E (^) - KP^Q) ■■ R < log(l + P/i),Q < Ph}, which 
becomes identical to the R-E region outer bound Cr_.e(P) (which is a box as defined earlier). 

Therefore, we know from the above discussions that only for the case of noise-free RF-band to baseband 
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processing (i.e., Case III), the PS scheme achieves the R-E region outer bound and is thus optimal. However, in 
practice, such a condition can never be met perfectly, and thus the R-E region outer bound Cr__e(-P) is in general 
still non-achievable with practical PS receivers. In the following, we will study further the achievable R-E region 
by the PS scheme for the more general case of MIMO channels. It is not difficult to show that if each receiving 
antenna satisfies the condition in Case III, the R-E region outer bound Cr_e(^) defined in ^ with G = H is 
achievable for the MIMO case by the PS scheme (with each receiving antenna to set p = 1). For a more practical 
purpose, we consider in the rest of this section the "worst" case performance of the PS scheme (i.e.. Case I in 
the above), when the noiseless antenna is assumed (which leads to the smallest R-E region for the SISO AWGN 
channel case). The obtained R-E region will thus provide the performance lower bound for the PS scheme with 
practical receiver circuits. In this case, since there is no antenna noise and the processing noise is added after 
the power splitting, it is equivalent to assume that the aggregated receiver noise power remains unchanged with 
a power splitter at each receiving antenna. Let pi with < pi < I denote the portion of power split to the EH 
receiver at the ith receiving antenna, I < i < N. The achievable R-E region for the PS scheme (in the worst case) 
is thus given by 

Cn-EiP) = U {iR^Q) -R< log 1^ + ^^J^HSH^K^J^IQ < tr{ApHSH^),tr{S) < P, S h 

0<Pi<l,Vi 

(13) 

where Ap = diag(pi, . . . , pn), and Ap = I - Ap. 

Note that the two points {Rma.x,0) and (0, Qmax) on the boundary of C^^g(P) can be simply achieved with 
Pi = 0,Vi, and pi = respectively, with the corresponding transmit covariance matrices given in Section HH 
(with G = H), similar to the TS case. All the other boundary points of C£^g(P) can be obtained as follows: Let 
H' = A^J'^H, G' = Ay^H, and n^^^{P,{pi}) denote the achievable R-E region with PS for a given set of 
Pi's. Then, we can obtain the boundary of 'R^_-^{P,{pi}) by solving similar problems like Problem (P3) (with 
H and G replaced by H' and G' , respectively). Finally, the boundary of C^-^{P) can be obtained by taking a 
union operation over different 7^^^(P, {pj})'s with all possible /jj's. 

In particular, we consider two special cases of the PS scheme: i) Uniform Power Splitting (UPS) with pi = p,\/i, 
and < p < 1; and ii) On-Off Power Splitting with pi G {0, l},Vi, i.e., pi taking the value of either or 1. For 
the case of on-off power splitting, let C {1, . . . , N} denote one subset of receiving antennas with pi = 1; then 
Q = {1, . . . , N} — Q denotes the other subset of receiving antennas with pi = 0. Clearly, Q and Q, specify the 
sets of receiving antennas switched to EH and ID receivers, respectively; thus, the on-off power splitting is also 
termed Antenna Switching (AS). 

Let 7^^^^(P, p) denote the achievable R-E region for the UPS scheme with any fixed p, and C^^|(P) be the R-E 
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region by taking the union of all Tl^^^{P, p)'s over < p < 1. Furthermore, let Tl^^{P, $7) denote the achievable 
R-E region for the AS (or On-Off Power Splitting) scheme with a given pair of Q, and Q. It is not difficult to see 
that for any P > 0, C^^KP) C CPS^(P), and n^^{P,n) C C^^_^{P),yn, while C^^l{P) = Cr^e(^) ^^r the 
case of MISO/SISO channel of H. Moreover, the following proposition shows that for the case of SIMO channel 
of H, C^lliP) = Cg^E(-P) is also true. 

Proposition 4.2: In the case of co-located EH and ID receivers with a SIMO channel H = h ^ C^^^, for any 
P > 0, CUP|(P) = CPS ^(P) = {{R,Q):R< log(l + m?P - Q)), <Q< \\hfP}. 

Proof: See Appendix IF] ■ 

D. Performance Comparison 

The following proposition summarizes the performance comparison between the TS and UPS schemes. 

Proposition 4.3: For the co-located EH and ID receivers, with any P > 0, C^%iP) ^ C^-e(P) ^ Cj^^e(^)' 
while CUP|(P) = Cl%{P) iff P < (1//12 - 1/hi). 

Proof: See Appendix IgI ■ 
From the above proposition, it follows that the TS scheme with the fixed power constraint performs worse than 
the UPS scheme in terms of achievable rate-energy pairs. However, the UPS scheme in general performs worse 
than the TS scheme under the flexible power constraint (without any peak power constraint), while they perform 
identically iff the condition P < (I//12 — is satisfied. This may occur when, e.g., P is sufficiently small 

(unlikely in our model since high SNR is of interest), or /12 = (i.e., H is MISO or SIMO). Note that the 
performance comparison between the TS scheme (with the flexible power constraint) and the PS scheme with 
arbitrary power splitting (instead of UPS) remains unknown theoretically. 

Next, for the purpose of exposition, we compare the rate-energy tradeoff for the case of co-located EH and ID 
receivers for a symmetric MIMO channel G = H = [I, 6;6, 1] with 9 = 0.5 for Fig. |7] and 9 = 0.8 for Fig. 
[H respectively. It is assumed that P = 100. The R-E region outer bound is obtained as Cr_e(P) with G = H 
according to Corollary 14.11 The two achievable R-E regions for the TS scheme with fixed vs. flexible power 
constraints are shown for comparison, and it is observed that C^_-^{P) C C^_^{P). The achievable R-E region 
for the TS scheme with the flexible power constraint P as well as the peak power constraint Ppcak = 2P is also 
shown, which is observed to lie between C^^g(P) and C^^g(P). Moreover, the achievable R-E region C^^|(P) 
for the UPS scheme is shown, whose boundary points constitute those of 7^^PS(P, p)'s with different p's from 
to 1. It is observed that C'^%{P) C C^P|(P) c C'^^^iP), which is in accordance with Proposition gj] Note 
that for this channel, the R-E region Cp^g(P) for the general PS scheme defined in ([T3] l only provides negligible 
rate-energy gains over Cj^P|(P) by the UPS scheme, and is thus not shown here. In addition, the achievable R-E 
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Fig. 7. Rate-energy tradeoff for a 2 x 2 MIMO broadcast system with co-located EH and ID receivers, and H = [1 0.5; 0.5 1]. 

region TZ^^{P,fl) for the AS scheme is shown, which is the same for Cl = {1} or {2} due to the symmetric 
channel setup. Furthermore, by comparing Figs. 0and|8] it is observed that the performance gap between Cp^^|(P) 
and Cj^^g(i-') is reduced when 9 increases from 0.5 to 0.8. This is because for this channel setup, hi = {1 + 9)'^ 
and /i2 = (1 - 9f, and as a result, the condition P < (I//12 - I//11) in Proposition |43] for C^%{P) = C^J^{P) 
will hold when —5- 1 for any P > 0. Finally, it is worth pointing out that in practical SWIPT systems with 
the co-located EH/ID receiver, under the high-SNR condition, the receiver typically operates at the "high-energy" 
regime in the achievable rate-energy regions shown in Figs. |7] and [H which corresponds to applying very large 
values of the time-switching coefficient a or the power-splitting coefficient p, i.e, a — )■ 1 and p — 1. 

V. Concluding Remarks 

This paper investigated the performance limits of emerging "wireless-powered" communication networks by 
means of opportunistic energy harvesting from ambient radio signals or dedicated wireless power transfer. Under a 
simplified three-node setup, our study revealed some fundamental tradeoffs in designing wireless MIMO systems for 
maximizing the efficiency of simultaneous information and energy transmission. Due to the space limitation, there 
are serval important issues unaddressed in this paper and left for our future work, some of which are highlighted 
as follows: 

• It will be interesting to extend the rate-energy region characterization to more general MIMO broadcast systems 
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8. Rate-energy tradeoff for a 2 x 2 MIMO broadcast system with co-located EH and ID receivers, and H = [1 0.8; 0.8 1]. 



with more than two receivers. Depending on whether the energy and information receivers are separated or 
co-located, and the broadcast information is for unicasting or multicasting, various new problems can be 
formulated for which the optimal solutions are challenging to obtain. 

For the case of co-located energy and information receivers, this paper shows a performance bound that 
in general cannot be achieved by practical receivers. Although this paper has shed some light on practical 
hardware designs to approach this limit (e.g., by the power splitting scheme), further research endeavor is still 
required to further reduce or close this gap, even for the SISO AWGN channel. 

In this paper, to simplify the analysis, it is assumed that the energy conversion efficiency at the energy receiver 
is independent of the instantaneous amplitude of the received radio signal, which is in general not true for 
practical RF energy harvesting circuits [21]. Thus, how to design the broadcast signal waveform, namely 
energy modulation, to maximize the efficiency of energy transfer to multiple receivers under practical energy 
conversion constraints is an open problem of high practical interests. 

Unlike the traditional view that the receiver noise and/or co-channel interference degrade the communication 
link reliability [22], they are however beneficial from the viewpoint of RF energy harvesting. Thus, there exist 
nontrivial tradeoffs in allocating communication resources to optimize the network interference levels for 
achieving maximal information vs. energy transfer. More studies to reveal such tradeoffs are worth pursuing. 
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Appendix A 
Proof of Proposition 12. II 

Without loss of generality, we can write the optimal solution to Problem (PI) in its eigenvalue decomposition 
form as Seh = VT,V^, where V G C^^^^^ VV^ = V^V = I, and E = diag{pi, . . . ,pm) with pi > 
P2 > • • • > Pm > and J2i=iPi — ^- Let G = GV = [^i, . . . tEIm]- Then, the objective function of Problem 
(PI) can be written as Q = tr{GSG^) = tr(GS;G^) = Y^fLiPiUiW^ < -PlI^ilP . where the equality holds if 
WqiW^ = maxj llgfjlP and pi = P,pi = 0,i = 2, . . . , M. Let V = [vi, . . . ,vm]- Since the (reduced) SVD of G 
is given by G = Uc^q'^Vg Section HIl we infer that WqiW^ is the maximum of all llgrjp's if and only if vi 
is the first column of Vg corresponding to the largest singular value of G, which is ^/gi. Hence, we obtain the 
optimal solution of Problem (PI) as Seh = Pviv^. The proof of Proposition 12. 11 is thus completed. 

Appendix B 
Proof of Theorem I3.1I 

The Lagrangian of (P3) can be written as 

L{S, A, /x) = log |/ + HSH^l + A (tr {GSG") - Q) - fi (tr(5) - P) . (14) 

Then, the Lagrange dual function of (P3) is defined as g{X, ji) = maxg'^Q L{S, A, /i), and the dual problem of (P3), 
denoted as (P3-D), is defined as min;^>o,^>o 9(A,/i). Since (P3) can be solved equivalently by solving (P3-D), 
in the following, we first maximize the Lagrangian to obtain the dual function with fixed A > and /U > 0, and 
then find the optimal dual solutions A* and ji* to minimize the dual function. The transmit covariance S* that 
maximizes the Lagrangian to obtain g{X*,f.i*) is thus the optimal primal solution of (P3). 

Consider first the problem of maximizing the Lagrangian over S with fixed A and fi. By discarding the constant 
terms associated with A and f.i in (fT4l ). this problem can be equivalently rewritten as 

maxlog 1/ + HSH^l - tr ((fil - XG^G) S) . (15) 

Recall that gi is the largest eigenvalue of the matrix G^G. We then have the following lemma. 

Lemma B.l: For the problem in ([TST i to have a bounded optimal value, ^ > Xgi must hold. 

Proof: We prove this lemma by contradiction. Suppose that fi < Xgi. Then, let S* = fivivf with /3 being 
any positive constant. Substituting S* into ([TSl l yields \og{l + fi\\Hvi\\'^) + fi{Xgi — p). Since H and G are either 
independent (in the case of separated receivers) or identical (in the case of co-located receivers), it is valid to 
assume that Hi/i^ilp > and thus the value of the above function or the optimal value of Problem ([TSl l becomes 
unbounded when /? — oo. Thus, the presumption that ^ < Xgi cannot be true, which completes the proof. ■ 
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TABLE I 

Algorithm for Solving Problem (P3). 

Initialize A > 0, > 0, p > Agi 
Repeat 

Compute S* using ( 117b with the given A and fj, 
Compute the subgradient of g{X,fi) 

Update A and /i using the ellipsoid method subject to > Xgi > 
Until A and jj, converge to the prescribed accuracy 
Set S* = S* 

Since Problem (P3) should have a bounded optimal value, it follows from the above lemma that the optimal 
primal and dual solutions of (P3) are obtained when fi > Xgi. Let A = fxl — XG^G. It then follows that A :^ 
with /i > Xgi, and thus A'^ exists. The problem in ([TSl l is then rewritten as 

maxlog|l + /i'5/i"^| -tr(A5). (16) 

Let the (reduced) SVD of the matrix HA'^/^ be given by iJA'^/^ ^ tjf^''^V^ , where tj G C^^'^^ V £ 
£MxT2^ f = diag(/ii, . . . , /itJ, with /ii > /i2 > . . . > /1T2 > 0. It has been shown in HTl under the CR setup 
that the optimal solution to Problem ([T6l ) with arbitrary A >~ has the following form: 

S* = A~^/^VAV^A-^/^ (17) 

where A = diag(pi, . . . ,pT:,), with pi = {I - l/hi)~^,i = l,...,T2. 

Next, we address how to solve the dual problem (P3-D) by minimizing the dual function g{X, ji) subject to 
A > 0, ;U > 0, and the new constraint ji > Xgi. This can be done by applying the subgradient-based method, e.g., 
the ellipsoid method ll20l . for which it can be shown (the proof is omitted for brevity) that the subgradient of 
g{X,fi) at point [X, fj] is given by [tr(GS*G^) — Q, -P — tr(S'*)], where S* is given in (fTTl). which is the optimal 
solution of Problem (ITSl ) for a given pair of A and fi. When the optimal dual solutions A* and fi* are obtained 
by the ellipsoid method, the corresponding optimal solution S* for Problem ([TS] ) converges to the primal optimal 
solution to Problem (P3), denoted by S*. The above procedures for solving (P3) are summarized in Table H The 
proof of Theorem 13.11 is thus completed. 

Appendix C 
Proof of Corollary 13. II 

Since H = , in Theorem O the (reduced) SVD of A'^/"^ with A = - X*G^G simplifies to 
^ ^ ^ ^ ^ ||A-i/2^f and vi = A-^/'^h/\\A-^/'^h\\. Thus, from © we have 
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S* = A-'/^vip.vi'A-'/^ (18) 
||A-V2/,|| V ||A-V2/,||2j ||A-V2/,|| ^ ^ 

Moreover, since r2 = 1 in this case, the maximum achievable rate is given by 

T2 T2 ^ + 

R* = J2^og{l + hip,) = J2{^ogChi)) = (21og(||A-i/2/i||)) . (21) 

i=l i=l 

From (Unil and dHJ, Corollary Othus follows. 

Appendix D 
Proof of Corollary 14. II 

Since G = H, from Theorem EH we have A = - X*G^G = - X*H^H y (i.e., fi* > X*hi). 
Recall that the (reduced) SVD of H is given by ii" = jy^r^V^^, with Th = diag(/ii, . . . , /it^), /^i > /12 > 
. . . > /it, > 0. Thus, it follows that A = /i*/ - X*H^H = Vh{i^*I - X*Th)V§, and A'^/^ = Vnif^*! - 
A*riy)-V2vg. Then, the (reduced) SVD of the matrix HA'^/^ is given by H{VHifJ-*I - A*r/^)Vg)-V2 = 
UHT]i'^{fJ-*I - A*rH)~^/Vg. Since in Theorem |331 the SVD of iJA^^/^ denoted by UT^^'^v" , we thus 
obtain U = Uh, f = Th{h*I - X*THy^, and V = Vh- From ©, it then follows that 

S* = A-^/^VAV^A-^/^ (22) 

= Vnif^*! - X*TH)-'/^V^VHAV^VH{fi*I - X*THr^^^VH (23) 

= Vnif^*! - X*Th)-^AV^ (24) 

^ Fi/EFg (25) 

where S = (^*/ - A*r//)^^A = diag(pi, . . . ,^^2)- Note that in Theorem |37T] A = diag(pi, . . . with 
= (1 - l/hi)+,i = 1,. . . ,T2, and f = diag(/ii,. . .,11x2) = ^nifJ-*! - A*r//)"^ Thus, we obtain that 

1 (^^^ 



/i* — X*hi \ h 
1 1 



i = i,...,r2. (27) 



^//* — X*hi hi 

Moreover, it is easy to verify that = TA. Since for Problem (P3), the maximum achievable rate is given by 



R* = ESi log(l + hiPi), it follows that 



R* = Y^log{l + h,p,). (28) 



1=1 



With (I25]l, and (l28]l, the proof of Corollary Ois thus completed. 
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Appendix E 
Proof of Proposition s. II 

Due to orthogonal transmissions for the EH and ID receivers in the TS scheme, we first show that the minimum 
transmission energy consumed to achieve any harvested power Q < Qmax in the EH time slot is equal to Q/hi 
regardless of a, as follows: From Section JI] (assuming G = H), it follows that the optimal S2 is in the form 
of qviv^ , where g > and vi is the eigenvector of the matrix H corresponding to its largest eigenvalue 
denoted by hi. To achieve Q, it follows from (ITTI ) that atr{HS2H^) = Q and thus q = Q/{hia). Thus, the 
minimum energy consumed to achieve Q in (fTTI ) is given by atr(S'2) = a ■ q = Q/hi, independent of a. 

With this result, in ([TTI ). the transmission rate R is given by (1 — a) log \I+HSiH^\ subject to {l — a)tr{Si) < 
{P — Q/hi). Due to the concavity of the log(-) function, it follows that R is maximized when a — 0, under which 
the optimal solution of Si can be obtained similarly as for Problem (P2). Thus, by changing the values of Q in 
the interval of < Q < Qmax and solving the above problem with a = 0, the corresponding maximum achievable 
rates as well as the boundary of Cj^^g(P) are obtained as given in (fT2l) for the case of flexible power constraint. 
Proposition 14.11 thus follows. 

Appendix F 
Proof of Proposition 14.21 

Since H = h = [hi, ... , h^]^ , for any set of p^'s with < < 1, the harvested power is equal to Q = 
P^f^i Pi\hi\^ . Clearly, Q <Q < \\h\\'^P. The equivalent SIMO channel for decoding information then becomes 
h = [^J\ — pihi, . . . , ^yl — p^hj^Y^ . Since for the SIMO channel, the transmit covariance matrix degrades to a 
scalar equal to P, the maximum achievable rate is given by (via applying the MRC beamforming at ID receiver): 

R = \og( I + \\hfp] (29) 



= \og{l + YSi-pi)\hi\''P^ (30) 

(N N \ 

1 + J]pi|/ii|2p (31) 
i=l i=l / 

= log{l + \\hfP-Q). (32) 

We thus have Cg^E(^) = {(^' Q) ■ R < log(l + {\\hfP - Q)), 0<Q < \\hfP}. Furthermore, since the above 
proof is valid for any /jj's and changing p from to 1 yields the value of Q = PpYl!i=i l^jP from to ||/i|pP, 
it thus follows that Cg!!|(P) = {(P,Q) : R < log(l + {\\hfP - Q)),0 <Q< P}, which is the same as 
C^^e(P). The proof of Proposition 14.21 is thus completed. 
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Appendix G 
Proof of Proposition |43] 

First, we prove the former part of Proposition O i.e., for any P > 0, C^%{P) C C^^|(P) C C^%{P). The 
proof of Cl%{P) C Cj^P|(P) is trivial, since the boundary of Cp^g(P) is simply a straight line connecting the 
two boundary points (0, Qmax) and (i^max, 0) (cf. Fig.|7]), and C^^|(P) is a convex set containing these two points. 
Next, we prove Cp^|(P) C Cj^^g(P),VP > 0, by showing that for any given harvested power < Q < Qmax, 
the corresponding boundary rate for C^^g(P), denoted by i?TS> is no smaller than that for denoted by 

i?UPS> i-e-> ^TS ^ PuPS> as follows: For any given Q, from the proof of Proposition 14.11 it follows that i?TS is 
obtained (with q = 0) by maximizing log \I + HSiH^\ subject to tr(S^i) < (P — Q/hi). On the other hand, for 
the UPS scheme, from the harvested power constraint ptr i^HSH^^ > Q, it follows that p > Q/{hiP) must hold. 
Note that i?uPS is obtained by maximizing log | / + (1 - p)HSH^\ subject to tr(5) < P. Let 5' = (1 - p)S. 
The above problem then becomes equivalent to maximizing \og\l + HS'H^\ subject to tr(S') < (1 - p)P. 
Since p > Q/{hiP), it follows that tr(S") < (1 - Q/{hiP))P = P- Q/hi. Thus, it follows that Pts > ^UPS- 
The former part of Proposition 14.31 is proved. 

Next, we show the latter part of Proposition |43] i.e., C^l%{P) = Cl%{P) iff P < (I//12 - ^/hi). Consider 
first the proof of the "if" part. For any < Q < (5max> since Pts = niax^'^ log|/ + HSiH^\ subject to 
tr(S'i) < (P — Q/hi) < (l//i2 — l/^i)> the optimal solution for this problem must be beamforming, i.e.. 
Si = {P — Q/hi)viVi with vi being the eigenvector of H corresponding to its largest eigenvalue hi, due 
to the WF power application given by (|3]l. Thus, it follows that Pts = log(l + hiP — Q). Consider now the UPS 
scheme. Suppose that S = Pvivf and p = Q/{hiP). It then follows that for UPS, the harvested power is equal to 
ptr [HSH^) = Q/{hiP) X {hiP) = Q, and the achievable rate Pups is equal to log | / + (1 - p)HSH^\ = 
log(l + (1 - Q/{hiP)) X hiP) = log(l + hiP -Q) = Pts- Thus, we prove that C^^_%{P) 5 C^MP)- Since 
from the proof of the former part of Proposition 14.31 we have that C^^|(P) C C^_^^{P), it thus follows that 
CUP|(P) = C^_\{P). The "if" part is proved. 

Second, we prove the "only if" part by contradiction. Suppose that Cg^|(P) = C^Ji^{P) for P = l//i2-l//ii+(5 
with arbitrary (5 > and thus P > (l//i2 — For any < Q < hi5, the corresponding Pts on the boundary 

of C^!.e(-P) is given by max^^ log |/ + HSiH^\, where tr(5i) = P - Q/hi > P - 5 = l//i2 - Since 
tr(S'i) > l//i2 — from the WF power allocation in <^ it follows that the optimal rank of S^i must be greater 
than one. Consider now the UPS scheme. Since CUP|(p) = C'^%{P), it follows that for the same Q as in the TS 
scheme, the maximum rate for the UPS scheme is Pups = Pts- From the proof of the former part of Proposition 
14.31 we know that to achieve Q with UPS, p > Q/{hiP) with the equality only when the transmit covariance is 
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rank-one. Furthermore, since 

i?UPS is equal to log |/ + HS'H^\ with tr(S") = {l-p)P < P-Q/hi, where the 
equality holds only if S' is rank-one. Since i?ups = ^TS> it thus follows that S' = Si and tr(S'') = P — Q/hi 
must hold at the same time. Since these two equalities require that S' have the rank greater than one and equal to 
one, respectively, they cannot hold at the same time. Thus, i?ups = ^TS cannot be true and the presumption that 
CUP|(p) = Cj^^"g(P) does not hold. The "only if" part is proved. 

Combining the proofs for both the "if" and "only if" parts, the latter part of Proposition 14.31 is thus proved. 
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